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Abstract 

The quark dynamics inside light mesons, except pseudoscalar ones, can be quite well described by 
a spinless Salpeter equation supplemented by a Cornell interaction (possibly partly vector, partly 
scalar). A mass formula for these mesons can then be obtained by computing analytical approxima- 
, tions of the eigenvalues of the equation. We show that such a formula can be derived by combining 

O 

■ the results of two methods: the dominantly orbital state description and the Bohr-Sommerfeld 

O 

D ■ quantization approach. The predictions of the mass formula are compared with accurate solutions 



of the spinless Salpeter equation computed with a Lagrange-mesh calculation method. 
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I. INTRODUCTION 



Semirelativistic potential models have been proved extremely successful for the descrip- 
tion of light mesons (mesons containing u, d, or s quarks). The main characteristics of the 
spectra of these mesons, except pseudoscalar ones, can be obtained with a spinless Salpeter 
equation supplemented with the Cornell interaction (a Coulomb-like potential plus a linear 
confinement) 

Numerous techniques have been developed in order to solve numerically with a great 
accuracy the semirelativistic equation. Nevertheless, it is always interesting to work with 
analytical results. Several attempts to obtain some mass formulae for hadrons were already 
performed. Some approaches rely on fundamental QCD properties jsl 0], but they are 
limited to the study of ground states of hadrons. In other works, the hadron masses are 
given as a function of some quantum numbers. They are based, for instance, on shifted 
large- iV expansion ( is the number of spatial dimensions) of the Schrodinger equation P| , 
a spectrum generating algebra j^, or a completely phenomenological point of view We 
will adopt here a different point of view by assuming that a semirelativistic potential model 
allows a good description of the main features of meson spectra. 

Recently, a new method to tackle this problem was developed: the dominantly orbital 
state (DOS) description, in which the orbitally excited states are obtained as a classical result 
while the radially excited states are treated semiclassically 0,0, 1^, [l^. A second method 
is the Bohr-Sommerfeld quantization (BSQ) approach, with which precise information can 
be obtained on the asymptotical behaviors of observables as a function of the quantum 



numbers 



12| |. We show here that a quite well accurate mass formula for light mesons. 



as a function of quantum numbers and parameters of a QCD inspired potential, can be 
obtained by combining the results of these two approaches. The idea is to calculate analytical 
approximate solutions of the equation assumed to govern the quark dynamics inside a meson. 
There is yet some uncertainties about the Lorentz structure of the interquark interaction. 
In this work, we will assume that the confinement potential is partly scalar and partly 
vector. A related work using a WKB approach was performed in Ref. but the Coulomb- 
like potential and the Lorentz structure of the confinement interaction were not taken into 
account. 

In order to test the validity of the mass formula, we compare its predictions with accurate 
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numerical solutions of the spinless Salpeter equation. These last ones are computed with a 



Lagrange-mesh calculation method 1^. This technique is here modified in order to handle 
semirelativistic equations with mixed scalar-vector potentials. 

In Sec. m the model Hamiltonian is presented with the two methods previously devel- 
oped to compute some analytical solutions. The mass formula is established in the case 
of symmetric and asymmetric mesons, with our without a constant term in the potential. 
In Sec. mil the mass formula is compared with accurate numerical solutions of the spinless 
Salpeter equation. Some concluding remarks are given in Sec. |TVl 

II. THE MODEL 

A. Model Hamiltonian 

Within the framework of a semirelativistic potential model, it is possible to describe the 
main characteristics of the spectra of light mesons . If the spinless Salpeter equation 
is chosen, instead of the Schrodinger equation, the quark-antiquark Hamiltonian is given by 
fwe use the natural units h = c = 1) 



H=^f + (mi + aiS{r)f + + {m^ + a2S{r)f + K(r), (1) 



where V{r) and S{r) are respectively the vector and scalar interactions |8|, and where p is 
the relative momentum between the quark and the antiquark. The vector p is the conjugate 
variable of the inter-distance r. As usual, we assume that the isospin symmetry is not 
broken, that is to say that the u and d quarks have the same mass (in the following, these 
two quarks will be named by the symbol n). The parameters ai and 0^2 indicate how the 
scalar potential S{r) is shared among the two masses mi and m2. A natural choice, used in 
this work, is to take 

^2 , mi 
ai = and 0:2 = . [2) 

mi + m2 mi + m2 

It is generally admitted that the short range part of the interquark potential is dominated 
by the one-gluon exchange process, which gives rise to a Coulomb term of vector type. The 
long range part is dominated by a confinement that lattice calculations predict linear in the 
interquark distance. As its Lorentz structure is not precisely known, we suppose here that 
the confinement is partly scalar and partly vector, as in Ref. P|. The importance of each 
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one is reflected through a parameter / whose value is for a pure vector, and 1 for a pure 
scalar. Consequently, the potentials considered here are given by 

Sir) = far (3) 

in which a is the usual string tension, whose value should be around 0.2 GeV^, and 

V{r) = il-f)ar-^ (4) 

in which k is proportional to the strong coupling constant a^. A reasonable value of k should 
be in the range 0.1 to 0.6. 

It is worth noting that it is not possible to describe the pseudoscalar mesons with a so 
simple potential. Spin contributions are very large in this sector as well as flavor mixing 
effects. An interaction stemming from instanton effects, which is not considered here, could 
exp.a,n tKe properties of tKese n,eson. QQ. Consequently, the psendoseala. n,esons cannot 
be described by our model. 

B. Semi-classical method 

Approximate analytical solutions of Hamiltonian (0) with potentials ©-0 can be ob- 
tained within the DOS approach. The idea of the model is to make a classical approximation 
by considering uniquely the classical circular orbits (lowest energy states with given total 
orbital angular momentum J), defined by r constant, and thus r = 0. The radial excitations, 
numbered with the quantum number u, are calculated by making a harmonic approxima- 
tion around the previous classical orbits. A detailed description of this method is given in 
Refs. Sayd. we just .eeaUhe. the n,a,n .suits. In the ease of a syn„netne .neson. 
nil = ^2 = ^1 the square meson mass is given by 

= aA{f)J + B{f)mVaJ + C{f)m^ + aD{f)K + aE{f){2u + l) + 0{J-^^^). (5) 

The coefficients A, B, C, D, and E are given by 

A{f) = j[t + 3(l-/)f , 

B{f) = y[(l + /)(3/-l)+t(l-/)], 

C{f) = -^[t(^ + f) + (l-/)(2/-l^ 



DU) = -|( + 3(1-/)| 



Eif) = m^J^, (6) 



where the auxihary functions s, t, and y are written 

s{f) = l-2f + 3f 



t{f) = ^s{f) + QP, 
g 

= s{f) + {i-f)t{fy 

The coefficients A, B, C, D, and E are monotonic functions of / and their ranges from 
/ = to / = 1 are 8 > A{f) > 4, < B{f) < 4v^, 8 > C{f) > 3, -4 < D{f) < -2^2, 
and 4-\/2 > > 4. Expression (jSj) is vahd for small values of m/ i/a and k, and/or large 
values of J. 

As this method relies basically on a classical approximation, it is not possible to calculate 
the zero point energy of the orbital motion. Thus, a mass formula cannot be obtained. 
Moreover, the dependence of the energy as a function of the radial quantum number u is 
calculated by making a harmonic approximation around classical orbits with high values 
of J. We cannot expect a good u behavior for small values of the angular momentum. A 
more serious flaw is that the method predicts a linear dependence of as a function of u 
whatever the form of the potential. So, we cannot be sure that the u dependence found is 
the more appropriate. A way to correct these drawbacks is to complete the previous analysis 
by a BSQ method. 



C. BSQ method 

The basic quantities in the BSQ approach are the action variables, 

Js=<l>Psdqs, (8) 

where s labels the degrees of freedom of the system, and where and Ps are the coordinates 
and conjugate momenta; the integral is performed over one cycle of the motion. The action 
variables are quantized according to the prescription 

J, = z/, + 1/2, (9) 
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where z/g (> 0) is an integer quantum number. This corresponds to a WKB expansion 



0Q,Q,Q 



hmited to the first order in h (see for instance 

The calculations for the angular momentum J, in the limit of great values for this quantum 

number, give simply the same J dependence for as in expression (jS]), but with J replaced 

by J+ 1/2. The calculations for the radial motion are more involved. A detailed description 

fl 

of the procedure is given in Ref. [I4I where the case of the Hamiltonian is studied for 
nil = ^2 and / = 0. We use here the same technique and we expand all expressions in 
powers of the meson mass M. Assuming that M'^ /a is large and J finite, and keeping only 
terms in M^, M, and 1/M, the integral (jH)) with Hamiltonian can be written, after 
tedious calculations. 




YZ X - 

with V the radial quantum number, and with 

X = M/ + 2m(l-/), 

Y = -4m^ + 2aK{l - /) + AamKf/M, 

Z = M(l - /) + 2m/ + aK{l - 2/)/M. (11) 

The expression above is valid only for f < 1/2. A similar equation exists for f > 1/2. It is 
now necessary to extract as a function of u in order to obtain an analytical result usable 
in a mass formula. If we assume that quantities m/^Ja and k are small, we can expand 
Eq. ()1()|1 in powers of these small parameters. The first order gives im/ ^fa = k = 0) 

^ aE\f){2v + 1) for > 1, (12) 

with 



1 2f f ^^Infii^^) for /<l/2, 

E'U) = 2. ] with if(/)= Vi-yp^; - (13) 

We have E'{Q) = 2%, E'{l/2) = 3n/2, and E'{1) = 4. This is in agreement with results 
obtained in Ref. for the case / = 0. 

The u square mass dependence obtained with this method is very similar to the one 
obtained with the DOS approach. But, the coefficients E{f) and E'{f) are different, as it 



is shown on Fig. ^ The approximations used to calculate these coefficients are also very 
different: E{f) is expected to give good results when J ^ u, while E'{f) is expected to give 
good results when z/ ^ J. 

By expanding the right hand side of Eq. (|Tn|l in powers of mj ^fa and k, we obtain, at 
the second order, 

'A(/) 



M^TT Mmf 1,^, A(f)f 



f 



1 



+ K. (14) 



where A(/) = 1 + y/1 — 2/ (note that this expression is well defined for / in the range [0,1]). 
We give these expressions for the sake of completeness. Nevertheless, as we can see below, 
the first order term is sufficient for our purpose. 



D. Mass formula 



Using results from the DOS approach and the BSQ method, we can write a square mass 
formula for light mesons composed of two identical quarks with a mass m as a function of 
the quantum numbers J and u, and the parameters of the potentials a, k, and / 

M2 = aA{f) (J+l/2) + 5(/) m^a{J+ l/2) + C(/) m^ + aD{f)K + a E'^'\f){2u + 1). (15) 

The coefficients A, B, C, D, and E are given by relations (jHI)-©, and the coefficient E' is 
given by relation (fT^ . 

Actually, this formula is an approximate expression for eigenvalues of the Hamiltonian (P). 
In principle, it is only valid for large values of J and/or large values of u. The coefficient 
E{f) will be preferred when J ^ u, while the coefficient E'{f) will be chosen when u ^ J. 
One can ask if the Eq. can give good meson mass estimations for realistic values of the 
potential parameters and for small values of the quantum numbers. In order to answer this 
question, and then to test the relevance of such a formula, it is necessary to compare the 
predictions of the formula with accurate numerical solutions of Hamiltonian (J]). This is the 
purpose of Sec. IIIII 



E. Addition of a constant term 



It is well known that it is necessary to take into account the contribution of a constant 
term in potential models in order to get the correct absolute values of the meson masses. In 
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principle it is possible to add a constant Ay to the vector potential and a constant A5 to 
the scalar potential. As in the case of the confinement potential, we can define a parameter 
g whose value is for a pure vector constant, and 1 for a pure scalar constant 

= (1 _ g)^ and A5 = gA. (16) 

A detailed discussion of the contributions of these two constants is given in Ref. pJi]. But it 
is very easy to see that the introduction a scalar constant is equivalent to a redefinition of the 
quark masses, while a vector constant simply shifts all meson masses. Finally, if M(m, A) is 
the mass of a meson containing two identical quarks interacting with the potential 0-0 
supplemented by a constant A given by Eq. |T6|l . then we have 

M(m,A) = M(m + ^A/2,0) + (1 -^)A. (17) 

Let us note that, as A is always negative in realistic potential models, we have m+gh./2 < m. 
Then, the formula ()17p gives a better approximation than without a constant, since it relies 
on an expansion for a small mass parameter. 



F. Asymmetric mesons 

The study of mesons containing two different quarks in the framework of the DOS ap- 



proach has been performed in Ref. [ll|. A formula for the square mass very similar to 
the expression has been found but with coefficients given only numerically. As these 
coefficients cannot be determined analytically, they are useless for a mass formula. 

Fortunately, one can verify experimentally that the mass of a ns meson is in very good 
approximation the arithmetic mean between the masses of the corresponding nn and ss 
mesons. Then, since the mass formula (fT3j) reproduces reasonably well the mass of a sym- 
metric meson, as we will see in Sec. IIIIl it is natural to try the following prescription 



X M(mi,mi) + M(m2,m2) . . 

M(mi,m2) = , (18) 

where M{m,m) is the mass of a symmetric meson given by Eq. and M(mi, 777,2) the 

mass of an asymmetric meson containing two quarks with masses mi and 7772. Equation ()18|l 
is then purely phenomenological. 

It is worth noting that the scalar potential and the scalar constant are equally shared 
between the two quarks in a symmetric meson. This is not the case for an asymmetric 
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meson. The relation ()18|) does not take into account this situation. We will verify below 
that this does not spoil the quality of the mass formula for asymmetric mesons. 

III. DISCUSSION OF THE MODEL 

In order to test the relevance of the square mass formula given by Eqs. (fT3j) . (|T7|). and 
()18|1 , we have compared its predictions with the exact eigenvalues of Hamiltonian (PJ) . They 
are obtained with a great accuracy by using a Lagrange-mesh calculation method for semirel- 
ativistic equations [l^. This method is described in the appendix, with the modification 
which must be made in order to handle the spinless Salpeter equation with mixed scalar- 
vector potentials. 

Realistic values for the quark masses and parameters of the potential are taken from a 
simple semirelativistic model described in Ref. m„ = 0.150 GeV, = 0.364 GeV, 
a = 0.203 GeV^, and k = 0.437 {vtinl \fa = 0.33 and m^j ^Ja = 0.81). In this paper, several 
constants are used, depending on the quark contents of the meson. Here we have just used 
the constant for nh meson, A = —0.599 GeV. The parameters are chosen to reproduce the 
main features of the meson spectra with f = g = 0. As, We have used these parameters 
in our square mass formula with varying values of / and g, we cannot obtain good spectra 
for all values of these two quantities. But the purpose of this work is simply to show the 
relevance of the mass formula. In Fig. El to Fig. we compare the exact eigenvalues of 
Hamiltonian with the predictions of the mass formula given by Eqs. fll5|) . (fT7|) . and ()18p. 

In Fig. 121 the square nn meson masses are plotted as a function of J and u for parameter 
values f = g = 0, the coefficient E{f) being used {E{f) and E'{f) are the most different for 
/ = 0). We can see that for small values of v or large values of J, the exact and approximate 
results are in good agreement. This situation is expected when the coefficient of the quantum 
number v is obtained with the DOS approach. For instance, for J = 0, the relative error on 
square mass increases regularly from 3.3% at z/ = 2 to 8.3% at = 9 (there are irregularities 
between i/ = and u = 2). In Fig. the same quantities are plotted, but the approximate 
results are calculated with the coefficient E'{f). In this case, for J = 0, the absolute error 
on square mass increases regularly from 0.120 GeV at z/ = to 0.570 GeV at z/ = 9, but in 
the same time, the relative error decreases from 24.5% to 3.1%. This is again the expected 
behavior for a coefficient of the quantum number u obtained with the ESQ method. In the 
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sector where the coefficients E{f) and E'{f) are different, it is interesting to choose the 
coefficient of the quantum number u following the values of the quantum numbers J and z/. 
With this constraint, the results of the square mass formula are quite good. 

When / = 1, the situation is simpler since E{f) and E'{f) are identical. In Fig. |3J the 
square nn meson masses are plotted as a function of J and u for parameter values / = 1 and 
g = 0. Again the agreement between exact and approximate results is good. The qualities 
of the mass formula seems to deteriorate with increasing values of u. Actually, for J = 0, 
the absolute error on square mass increases regularly from 0.148 GeV at z/ = 2 to 0.674 GeV 
at u = 9, but in the same time, the relative error decreases slowly from 5.9% to 5.5% (there 
are irregularities between 1^ = and u = 2). 

Finally, we have tested the formula (fTH|) by calculating the square ns meson masses as 
a function of J and u for parameter values / = 0.5 and g = 0.5 with the coefficient E'{f). 
Results are plotted in Fig. We can see that the prescription works very well. 



IV. CONCLUDING REMARKS 



The main features of the spectra of light mesons, except pseudoscalar ones, can be re- 
produced with a spinless Salpeter equation supplemented with the Cornell interaction 
We have shown that the eigenvalues of this simple Hamiltonian can be obtained, within a 
few percents of relative error, by a mass formula. This relation gives the square mass of 
light mesons as a function of quantum numbers J and u and the parameters of the po- 
tential. One could expect that the formula is only valid for high values of these quantum 
numbers and/or for very small values of parameters m/y/a and k. Outside these limits, we 
have remarked that the formula only differs slightly from exact results for physical values of 
parameters. This allows to apply our approximation for physical situations and to compare 
our calculations with experiment. 

The values of the string tension a and of the strength k, of the Coulomb-like potential can, 
in principle, be computed from lattice calculations. The values of the constituent masses and 
of the constant potential are more difficult to obtain, as well as the quantities / and g. These 
parameters, in particular the scalar- vector mixture in the confinement, could be determined 
by a fit of our square mass formula on available experimental data. Unfortunately, some 
uncertainties on experimental meson masses, leading to large intervals of possible values for 
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the parameters, make difficult such a determination for the moment |10|. 

The meson mass formula obtained here shares some similarities with other formulae, but 
,t .Ue. on d«e..e,. bas,s: In.ead of .eMng o. a spec... ,en.at,„g a.,eb.a Q o. on 
a completely phenomenological point of view it is assumed here that a semirelativistic 
potential model allows a good description of the main features of meson spectra, as it is the 
case in Ref. jl3|. The meson mass dependence on some usual parameters of semirelativistic 
potential model is then obtained with a good approximation. 

In our study, we completely neglect the quark spin. In usual models the spin-dependent 
part of the potential is the hyperfine interaction stemming from the_pne-gluon exchange 
interaction (and may be from the vector part of the confinement) |2l|. In other models, 
this spin-dependent part stems from an instanton induced interaction jlll. Il5|. Except for 
pseudoscalar mesons, in both cases, the contribution of the spin to the meson masses is small 
or vanishing with respect to orbital or vibrational excitations. Then, despite the absence 
of spin dependence in our Hamiltonian, our mass formula can describe a large sample of 
mesons. 

NUMERICAL METHOD 

The semirelativistic Lagrange-mesh method can be used to compute the eigenvalues and 
the eigenf unctions of the following spinless Salpeter Hamiltonian 

if = y^^Tmf +K(r). (19) 



A detailed presentation of the technique is given in Ref. [1^. Only the main points of the 
method are given here. A variational calculation is performed with the trial state 

l^) = E^.I/.) where (f |/,.) = (20) 

The coefficients Cj are linear variational parameters and the scale factor h is a. non-linear 
parameter aimed at adjusting the mesh to the domain of physical interest. The functions 
fj{x) are such that 

Mx,) = x;'/'6,,. (21) 
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The numbers Xi, which are the zeros of a Laguerre polynomial of degree N, and the numbers 
Aj are connected with a Gauss quadrature formula 

^0 k=l 

At the Gauss approximation, {fi\fj) ~ Sij, the potential matrix elements are simply given 

by 

{mir)\f^)^Vihx,)S,,. (23) 



The computation of the matrix elements {fi\\/p'^ + rn'^lfj) is obtained from the calculation 
of the matrix elements + m'^\fj). These last quantities can be easily obtained at 

the Gauss approximation with analytical formulae Q|- Despite the use of an approximate 
quadrature rule, a very high accuracy can be attained with a small number of basis states. 
When a Hamiltonian of the form is considered, the matrix elements of the operators 



P + (^fc + 0!kS{r)) must be computed. Again, they can be obtained from the matrix 
elements of the square of the operators. From Eq. (j^^ . we have 



{Mf + {rrik + akS{r)f |/,) ^ {f,\p' + ml\f,) + {2mkakS{hx,) + alS{hx,f) 6,,. (24) 

We have checked that this procedure allows the computation of the eigenvalues and the 
eigenfunctions of the spinless Salpeter Hamiltonian (Q) with a high accuracy. 
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FIG. 2: Square masses of nn mesons as a function of u and J quantum numbers, with param- 
eters of Ref. Q], and for g = and / = 0. Solid lines join exact solutions of the spinless Salpeter 
equation. Approximate results from Eq. H15|) with the coefficient E{f) are indicated by circles. 
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FIG. 3: Same as Fig. [2 but for Eq. ((TSJ) with the coefficient E'{f). 




FIG. 4: Square masses of nn mesons as a function of u and J quantum numbers, with param- 
eters of Ref. ^1, and for g = and / = 1. Solid hues join exact solutions of the spinless Salpeter 
equation. Approximate results from Eq. (|15j) are indicated by circles {E'(l) = E{1)). 
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FIG. 5: Square masses of ns mesons as a function of and J quantum numbers, with parame- 
ters of Ref. Q], and for g = 0.5 and / = 0.5. Solid lines join exact solutions of the spinless Salpeter 
equation. Approximate results from Eq. H18() with the coefficient E'{f) are indicated by circles. 
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